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Abstract 

Using a local effective potential to account for 
electron exchange, we have analyzed the zero-energy scat- 
tering of ortho-positronium by helium atoms. The effective 
potential was chosen so that the results obtained in the 
static approximation agreed with the static -exchange results 
of Fraser and Kraidy. Correlation was introduced in two 
different approximations which agree well. We find a de- 
crease in the scattering length a from I.72 a Q to 1.39 a Q j 
and an increase of more than a factor of 2 in the quenching 
rate of the positronium (from Z eff = 0.0^2 to 0.10), as 
compared with the static exchange approximation. This 
gives improved agreement with the 


Z eff 


0.18. 


experimental value 


1. Introduction 


The system consisting of one atom of helium and one 
of the long-lived ortho-positronium is accessible to ex- 
periment through observation of the density dependence of 
the annihilation rate of positrons in helium gas. The 
positron can be "picked off" or quenched during collision 
with a helium atom and this quenching rate is therefor^ 
proportional to the helium density. It measures the joint 
probability that the positron be at the location of one of 
the helium electrons and also form a singlet spin state 
with that electron. This pickoff probability is convention- 
ally presented in terms of the effective number of electrons 
in the singlet state seen by the positron, Z e ^. This 
quantity is, in general, energy dependent, and is obtained 
from the scattering wave function t 

More specifically, if Y(x, £l, rp, rp) is the Ps-He 
scattering function, x being the coordinates of the posi- 
tron and ri, Tz, 1*3 the coordinates of the electrons then 

•■v* 

one can write 


X J $ (ji: “ ) 


(l) 



where is the projection of the wave function Y on to the 
singlet spin function involving the positron and the i th 
electron. 

Fraser and Kraidy (1966) and Fraser (1968) have calcu- 
lated Z e;ff and. the scattering length a in the static 
exchange approximation. They assumed a trial function of 
the form 

* (^1.^3 ~ fixes') } ^ 

where ^ = £(x + r^), = r^- x , 0 and $ are the posi- 

tronium and helium ground state wave functions respectively. 

is an operator which antisymmetrizes the wave function 
with respect to the three electrons, and a, j3 represent up 
and down spin respectively. They del :.rmined X "by solving 
numerically a certain integro -differential equation obtained 
variationally, and with the asymptotic normalization 
X(R) ~ exp (ik.R), they used (l) to compute Z eff> . For zero 
energy collisions their results were = 0.042, and the 

scattering length a * 1.72 a Q .* This value for Z eff is in 


3 


*Qhese numbers are slightly different from those published 
by Fraser and Kraidy (1966). The change is due to a dif- 
ferent method of treating the target wave function (Kraidy, 
1969, private communication). 

strong disagreement with the best experimental value of 
0.l80 4 0.016, deduced by Fraser (1968) from the experi- 
mental data of Beers and Hughes (1968). Earlier experimen- 
tal results for Z eff . range from 0.118 4 0.011 (Duff and 
Heymann, 1962) to 0.25 ± 25$ (Roellig and Kelly, 1967). 

This static exchange approximation satisfies the con- 
ditions on scattering bounds (Spruch and Rosenberg, i960) 
and hence the scattering length obtained is larger than the 
true value. Although no rigorous bound theorems for Z eff 
are known, it is roughly correct that a large, positive 
scattering length must correspond to a small Z^^. One 
expects, then, that any improvement in the trial wave func- 
tion (2) leading to a decrease in a will also tend to 
increase Z eff* 

Barker and Bransden (1968, 1969) examined the effect of 

van der Waals’ long range forces on Z^^. Such forces are 

attractive, reduce a and are expected to increase Z 

“ eff 

They found, however, that an increase of only about 30$ in 


k 




Z __ results from this cause, 
ell 

In this paper we take the view that, since the strongest 
force acting between the two atoms is the exchange force be- 
tween the closed electron shell in helium and the electron 
of positronium, short-range correlation induced by this 
force will be the dominant correction to the trial function 
(2). Since Ps is much more polarizable than He, most of 
this correlation will affect the former, while the latter will be 
almost undistorted. In the next section, we outline a 
model embodying these considerations which, although simple, 
seems realistic enough to describe the system fairly well. 

2. Description of the model 

Consider the following model: assume the Ps electron 

to be distinguishable from the He electrons, and introduce 
an effective local interaction between the helium atom and 
this electron to represent the effect of the exclusion 
principle. Specifically we replace the exact Hamiltonian 
of the system (in atomic units, with energies in Rydbergs) 

>y the following: 

H • H. - 5 P,‘ + V.t" 1 ’" (» 
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where H q = - 2V p 2 - | + H(l,2) 

where coordinates without indices refer to electron 3, now 
treated as distinguishable, and H(l,2) is the Hamiltonian 
for the helium atom alone. The ground states, 0(p) of Ps, 
and t( r i> r 2) of He, satisfy the equation 

(H-e.) f(f) O = o y M 

and the model potential depends on two parameters, V Q and a. 
We have omitted the electrostatic interaction with the 
helium charge density and, since the helium atom is not 
distorted in this model, no van derWaals* force occurs 
either; these two should have small effects which would 
partially compensate. 

To fix the values of V Q and a appearing in(3) we use 
our model without distortion of the positronium, and demand 
that the results agree with Fraser and Kraidy (1966) and 
Kraidy (private communication) at zero energy. That is we 
assume a wave function analogous to (2): 

f -- %(*) f>(r) HW) { ' J) 
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where we have suppressed the spin indices. Using the 
Hamiltonian of (3) in a standard variational principle and 
recalling (4) we obtain the equation satisfied by x(R) in 
the static model approximation: 


[- 1 V * *u i (*)]k(«') 


* o 


(6) 


where 




(7) 


since 


r-- 


Equation (6) was integrated numerically to give both x(R) 
and the scattering length a. 

Using (j> (f) = (-{[) (8) 

and the approximate helium wave function 

, 2 * 71 , (9> 
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we compute the effective electron nunber Z from the 

eff 

expression 


z eff * 2 do) 

The parameters V Q and a were adjusted to give the values of 
a and obtained by Fraser and Kraidy. We found 
V Q » 39 and a = 2.9 gave sufficiently good agreement. 

5. Consequences of the model 

With the parameters of the model defined, we will next 
calculate, by two different techniques, the values of a 
and Z eff which result from the Hamiltonian of (3). We will 
consider the zero- energy case only. 

3.1 Differential Equations 

To the trial function of (3) we add a closed ch anne l 
term giving 

* [x&O *(f) (11) 
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whore G is a correlation function to be chosen later* By 
a method like that leading to equation ( 5 ), we now obtain 
the pair of coupled differential equations 


U?; + ‘Ufa)? * ° 

[.i n7 ;*wu„£ *h S 4.5Q f m1 

4- % = O • 


If we define the notation 

<0 * Hr) { Hr) (13) 


and assume, without loss of generality, that < G ) * 0, 
we can write the quantities appearing in ( 12 ) as follows: 


ft,- <» , «,• , «j ■ 4 .? v >- 


n * 0 "> , w* 9 )> , 
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-ar 

where V = V Q e is the effective potential appearing in 
(3). This notation is based on that used in two previous 
papers (Drachman 1968 , 1969 ) • Once the solution of equation 
(12) is obtained, Z eff . can be calculated by replacing X in 
(10) by (x + FG). 

Guided by perturbation theory in a closure approxi- 
mation, we have used an exponential form for the correlation 
function. If 

y I? 

•m 

- 

which satisfies the condition (g) = 0, and depends on the 

parameter y . In figure 1 we show the dependence of a 

and Z „ on y . The best (minimum) value of a occurs at 
err — 

y = 0.824 , where a = 1.42 a and Z __ = 0.100 . 

' — o eff 

3.2 Kohn Variational Principle 

In this model scattering problem, Kohn’s variational 
principle may be written as 

a $ (15> 


S+. - 4*. [- 


we take 


9 * * 


10 


where H is given hy (3) and the trial function has the 
asymptotic form 


f T ~ t(r) H *** 0 <*> 

as R - 00 . Note the factor 2 outside the integral in (15). 
Ibis does not appear in Kohn's variational principle for 
positron-atom scattering and arises from the use of the co- 
ordinates R and P in place of x and r . As already mention- 
ed, this system satisfies the conditions for the variational 
approximation a y to be an upper bound to the exact scattering 
length (Spruch and Rosenberg, i960) . Using the separable 
form of trial function Y given by (5) with 


#(*)• -*] 


where a^, is the trial scattering length and 6 is a non- 
linear variational parameter, we found that a^ has the 
minimum value 1.763 a Q when 6 = 0.6 . We then modified 
X(R) by adding on terms 



JL 




z 
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until a converged. This was with N = 4 and then a was 
v v 

1.722 a Q which is in agreement with the value obtained by 
numerical integration in the previous section and thus with 
Fraser and Kraidy. Furthermore when Z eff , is calculated 
using (10), with the parameters and a^ determined by the 
variational calculation, we obtain the value Z = 0.042, 
thus verifying that our model reproduces the static exchange 
results of Fraser and Kraidy. 

We now introduce correlation into the calculation by 
replacing x(R) by a function X'(R,P,r) thus: 

r)* rfiTfi-e**) -*} 

A.»t I 

where 3, y are additional non-linear parameters and 

i 0 . This is a generalized Hylleraas type ex- 
pansion in the coordinates R, P and r = |r + % p | and is 
anologous to using the coordinates (ri, r 2 , ri 2 ) in, for 
example, a positron-hydrogen atom scattering problem. 

(See, for example, Schwartz, I 96 I. ) This is an effective 
way of including the distortion of the Ps atom. 


12 



The number of terms N was chosen to include all terms 


such that 

+ m. + n. 
i i 1 

where u> took the values 1 to 6. This corresponded to N 
taking the values 10, 20, 35* 56, 84. At this stage of 
the calculation all the integrals were carried out analyt- 
ically. 

With N = 35* the three non-linear parameters 6, 0, y, 

were varied until the minimum value of I.389 a Q was obtained 

for a v . This was with 6 = 0.7, 0 = -0.2, V = 0.8. This is 

close to the value of 1.42 a obtained by the differential 

o 

equation method, but is lower because a more general cor- 
relation function was used. Increasing the number of terms 
K to 84 did not reduce a y below this value, and further- 
more, with N = 84, the value for was Insensitive to small 
changes in the non-linear parameters. We conclude that 
I.389 a Q is the least upper bound obtainable for the scatter- 
ing length with this type of trial function. Table I shows 
the convergence of a^ as N is increased. 

We now substituted this variationally determined 35- 
term wave function into equation (10) and calculated 
We now used numerical integration and the value obtained 
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for Z eff was O.O98, again in good agreement with the result 
of Section 3.1. 

We could only have calculated the Z eff integrals analyt- 
ically by taking y » 0 and restricting n^ to be even. In this 
case we found that did not converge even after 80 terms 
(corresponding to u> a 7) and was still as large as 1.463 a Q . 

We abandoned this line of enquiry. However, this point has an 
interesting analogy in positron-hydrogen atom scattering. In 
the present coordinate system, (R, p, r), r corresponds to the 
"interelectronic" coordinate r i2 used in Hy Her aas -type wave 
functions for positron-hydrogen scattering and it is well known 
that an expansion which includes all powers of ri 2 converges 
much faster than an expansion in Legendre polynomials of the 
cosine of the angle between £1 and r^. This polynomial ex- 
pansion essentially restricts the power of ri 2 to even values. 
We have seen this demonstrated again in the present problem. 

4. Conclusions 

Within the restricted framework of our simple model we 
have examined the zero-energy scattering of ortho-positronium 
by helium. Short-range correlation terms representing dis- 
tortion of the Ps atom were effective in bringing about a 
large Improvement in the pickoff annihilation rate, sub 



compared with experiment. At the same time, the scattering 
length was significantly reduced. 

A model of this kind serves only to indicate a direc- 
tion for future work. In this case, it suggests that a 
variational calculation with short-range correlation terms 
should prove successful in treating the fundamental Ps-He 
problem, as an extension of the earlier work of Fraser 
and Kraidy (l 966). 

We emphasize the importance of the Ps-He system: it 

is the next step in complexity after the e + -He system, and 
is both experimentally and theoretically tractable. 
Measurements of the quenching rate are already available, 
as noted above, and the interesting experiments on cavities 
in liquid helium (Roellig and Kelly, 1967) offer the pos- 
sibility of measuring the scattering length. We hope that 
the results reported here will encourage others to examine 
both theory and experiment more fully. 
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Table 1. 


a and Z __ for different values of N 
v err 


N 

4 

10 

20 

35 

56 

84 

a 1.448 

V 

1.396 

1.391 

1.389 

1.389 

1.389 

Z eff 

0.090 

0.096 

O.O96 

0.098 
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Caption for Figure 1 


a and Z e ^ as a function of y f calculated by the 
differential equations method. 
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